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ABSTRACT 

We use high signal-to-noise spectroscopy for a sample of 232 quiescent galaxies in the Shap- 
ley Supercluster, to investigate how their stellar populations depend on velocity dispersion 
£Nj ■ {cf), luminosity and stellar mass. The sample spans a large range in velocity dispersion (30— 

300 km s -1 ) and in luminosity (Mr from -18.7 to -23.2). Estimates of age, total metallicity 
(Z/H) and a-element abundance ratio (a/Fe) were derived from absorption-line analysis, using 
single-burst models of Thomas and collaborators. Using the Rose Ca II index, we conclude 
, that recent star-formation ("frosting") events are not responsible for the intermediate ages 

observed in some of the galaxies. Age, Z/H and a/Fe are correlated positively with velocity 
dispersion, but we also find significant residual trends with luminosity: at given a, the brighter 
galaxies are younger, less a-enriched, and have higher Z/H. At face value, these results might 
suggest that the stellar populations depend on stellar mass as well as on velocity dispersion. 
However, we show that the observed trends can be reproduced by models in which the stellar 
populations depend systematically only on a, and are independent of stellar mass M*. For 
age, the observed luminosity correlation arises because young galaxies are brighter, at fixed 
M* . For metallicity, the observed luminosity dependence arises because metal-rich galaxies, 
at fixed mass, tend also to be younger, and hence brighter. We find a good match to the ob- 
t> ' served luminosity correlations with Age oc <j +oao , Z/H oc cr +0 35 , a/Fe oc cr +a2 °, where the 

, slopes are close to those found when fitting traditional scaling relations. We conclude that the 

star formation and enrichment histories of galaxies are determined primarily by the depth of 
their gravitational potential wells. The observed residual correlations with luminosity do not 
imply a corresponding dependence on stellar mass. 

OO 
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! '~j 1 INTRODUCTION average harbour older stars, with higher total metallicity, and with 

r> higher abundance of a elements relative to Fe, as compared to less 

I_h 1 Numerous studies have explored how the stellar content of passive . . . 

J3 ■ , . JT1 , , • „ » »i , massive galaxies, 
galaxiesrj depends systematically on mass . Along the passive se- 
quence, more massive galaxies are observed to be redder, and have While the strong mass-dependence of passive galaxy proper- 
stronger metal absorption and weaker hydrogen features in their ties is widely appreciated, it is less clear which aspect of "mass" 
spectra (e.g. Faber 1973; Bower, Lucey & Ellis 1992). Physical ex- is really dominant in establishing the observed correlations. From 
planations for these trends have been sought using increasingly so- simple physical arguments we might expect metallicity to depend 
phisticated stellar population synthesis models (e.g. Worthey 1994; on the retention of gas within the galaxy halo during a galactic wind 
Thomas, Maraston & Bender 2003; Schiavon 2007), and ever im- phase (e.g. Arimoto & Yoshii 1987), and to be primarily correlated 
proving observational datasets (e.g. Trager et al. 2000; Kuntschner with escape velocity, and thus with velocity dispersion a. Alterna- 
et al. 2001; Caldwell et al. 2003; Nelan et al. 2005; Thomas et al. lively, if the amount of metals and stars produced both reflect the 
2005; Bernardi et al. 2006; Graves et al. 2007). These studies have overall efficiency of star-formation, then metallicity might instead 
converged on a scenario in which more massive passive galaxies on be correlated primarily with stellar mass, and hence with luminos- 
ity. Because luminosity and velocity dispersion are mutually cor- 
related through the Faber-Jackson (1976, FJ) relation, and because 
1 By the terms passive/quiescent galaxies, we loosely include many stellar mass is inferred rather than measured directly, it is not triv- 
partially-overlapping definitions, i.e. galaxies which have ceased forming ial to distinguish the relative importance of these mass proxies in 
stars, lie in the red envelope of the optical colour distribution, have little or driving the relations. By exploiting the substantial scatter in the FJ 
no line emission in their spectra, have early-type morphology, etc. relation, a number of previous works have indicated that luminosity 
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and velocity dispersion are not interchangeable as predictors of the 
stellar populations (e.g. Bernardi et al. 2005; Gallazzi et al. 2006). 
In particular, Graves, Faber & Schiavon (2009) recently analysed 
a sample of galaxies from the Sloan Digital Sky Survey (SDSS), 
concluding that important correlations with stellar mass persist af- 
ter accounting for the trends with velocity dispersion. 

In this paper, we address these issues using results from high 
signal-to-noise spectroscopy for galaxies in the Shapley Superclus- 
ter. The observations were described in detail by Smith, Lucey & 
Hudson (2007, hereafter Paper I), which also presented measure- 
ments of absorption line strengths for a sample of quiescent su- 
percluster member galaxies. Here we employ the data to deter- 
mine ages, total metallicities (Z/H) and a-element abundance ratios 
(a/Fe) for each galaxy in the sample, and analyse their dependence 
on luminosity and velocity dispersion. 

Section [2] provides a brief review of the sample definition 
and spectroscopic data, and describes the estimation of single- 
burst equivalent stellar population parameters. We also provide 
here a test for consistency with the broadband colours, discuss 
the evidence against widespread recent secondary bursts of star- 
formation, and describe the corrections for metallicity gradients. In 
Section [3] we present the correlations of age, Z/H and a/Fe with 
velocity dispersion and luminosity. In Section [4] motivated by the 
significant correlations of the stellar population parameters with lu- 
minosity, at fixed velocity dispersion, we investigate whether these 
require an underlying correlation with stellar mass. Section[5]com- 
pares the results to previous findings, especially focusing on the 
work of Graves et al. (2009), and our main conclusions are summa- 
rized in Section[6] 

The stellar population parameters as determined here have 
been used by Rawle et al. (2008a) to study the dependence of 
ultra-violet colours on the stellar populations, and by Gargiulo et 
al. (2009) to analyse stellar population effects in the Fundamen- 
tal Plane residuals. An alternative analysis of a subset of the Paper 
I data, using the method of Graves & Schiavon (2008) has been 
presented by Smith et al. (2009a). A more thorough exploration of 
complex star-formation history models is provided by Allanson et 
al. (2009). A later paper will treat the correlations between stellar 
population parameters, e.g. the age-metallicity-mass relation. 

Throughout this paper, we adopt cosmological parameters 
(QM,i^A,h) — (0.3,0.7,0.7). For reference, at the redshift of 
Shapley (z — 0.048), one arcsecond corresponds to 0.97 kpc, and 
the distance modulus is m — M = 36.65. 



2 STELLAR POPULATION PARAMETERS 
2.1 Data 

Paper I provides a full description of the observations and data pro- 
cessing. Here, we provide only a summary of the main points. 

The initial galaxy sample was drawn from photometric cata- 
logues of the NO AO Fundamental Plane Survey (NFPS, Smith et al 
2004), covering the central 40x40 arcmin 2 region in each of three 
clusters, Abell 3556, Abell 3558 and Abell 3562, in the core of 
the Shapley Supercluster. Spectra were obtained using the fibre-fed 
dual-beam AAOmega spectrograph at the Anglo-Australian Tele- 
scope. The fibres sample an aperture of 2arcsec diameter. Within 
the region covered by the NFPS catalogue, we observed some 
60 per cent of all R < 18 galaxies, with little dependence on mag- 
nitude. In the blue arm of the spectrograph, the spectra cover the 
main Lick indices, up to Fe5406 at the red end, with a spectral res- 



olution of 3.2 A FWHM. The red arm was used to record the Ha re- 
gion, at 1 .9 A resolution, for the purpose of detecting nebular emis- 
sion. The total integration time was ~8 hours per galaxy, result- 
ing in high signal-to-noise ratio (median S/N~60A _1 at 4400- 
5400 A). 

The primary parameters measured from the spectra are the 
redshift, cz, the velocity dispersion, a, equivalent widths of various 
emission lines after removing the stellar continuum, and the Lick 
absorption line indices. For some galaxies, the velocity dispersion 
is "unresolved", i.e. indistinguishable from zero in our spectra. The 
absorption line indices were transformed to the spectral resolution 
of the Lick system and corrected for the effects of velocity broad- 
ening. A final sample for stellar population analysis was defined 
which excludes galaxies with Ha emission equivalent width above 
0.5 A (which indicates likely contamination of the Lick Balmer-line 
indices), and those which lie outside of the adopted redshift range 
for the supercluster (cz — 11670 — 17233 km s" 1 ). The final sam- 
ple comprises 232 galaxies, of which 198 have measured velocity 
dispersion and 34 are unresolved. The analysis presented in this pa- 
per is based on the line-strengths and supporting data provided in 
tables 2^1 of Paper I. 

Note that the sample was not selected explicitly according to 
morphology at any stage, and is therefore representative of the pas- 
sive population overall. Although no selection on colour was ap- 
plied, our Ha selection very effectively restricts the sample to the 
red sequence in the B — R colour magnitude relation (see figure 1 
of Paper I). 

2.2 Linestrength diagrams 

This and the following section describe the process of transform- 
ing the set of index measurements for each galaxy into a set of 
physical parameters, by comparison with stellar population mod- 
els. This "grid-inversion" approach has of course been applied by 
many other studies in the field, including some fairly sophisticated 
implementations (e.g. Proctor, Forbes & Beasley 2004; Kelson et 
al. 2006; Graves & Schiavon 2008). Here, the absorption line data 
will be interpreted using simple stellar population (SSP) models 
by Thomas and collaborators (Thomas, Maraston & Bender 2003; 
Thomas, Maraston & Korn 2004, hereafter collectively TMBK). 
Constraints on non-SSP formation histories will be discussed in 
Section|231 

For the TMBK models, just three non-redundant index mea- 
surements are sufficient to invert the model grid and recover esti- 
mates of log(Age), [Z/H] and [a/Fe]. If more indices are available, 
a \ 2 minimization provides an estimate of the SSP which best re- 
produces the observations. Such redundancy may be desirable both 
to suppress random errors, and to average over systematic error 
sources (e.g. the dependence of indices on chemical abundances 
beyond those allowed to vary in the models). For age-sensitivity, at 
least one index should measure one of the Balmer lines. In princi- 
ple, any two metal lines with different a-sensitivity could be used 
to distinguish Z/H from a/Fe effects. With the advent of improved 
models which follow the light-element abundances separately (e.g. 
Schiavon 2007), it is important to be clear which elements are con- 
tributing to the measured a/Fe. Mgb5177 is the only index em- 
ployed as an a-abundance indicator, and although we will use the 
notation [a/Fe] in the context of the TMBK models, our measure- 
ments primarily reflect [Mg/Fe] . 

Our choice of Fe-sensitive index is limited because the widely- 
used Fe5270 and Fe5335 features are contaminated by 5577 A 
night-sky emission line in most of our spectra. The other red Fe 
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Figure 1. Index-index planes showing the a/Fe-sensitive index pair Fe5015 and Mgb5177. Each panel shows galaxies within a given interval in velocity 
dispersion. For clarity, we omit the ~4 per cent of galaxies having S/N < 20, and two other galaxies with large error bars in Fe5015 due to incomplete data. 
The model grids are from Thomas et al. (2003, 2004), with zero-points re-calibrated from the most massive galaxies (see Section |2\3"V The grids show lines 
of constant [a/Fe] = 0.0, +0.2, +0.3, +0.5 (red) and of constant [Z/H] = —0.33, 0.00, +0.35, +0.67 (green). For this comparison, the grids are drawn for a 
fiducial age in each panel, as indicated in the legend. The fiducial ages are estimated from visual examination of Figure[2] 



index, Fe5406, may be unreliable because it is close to the spec- 
trograph dichroic cut-off. Among the bluer lines, Fe4668 is also 
known as C24668 because it is dominated by the Swan C2 band. 
Since the abundance behaviour of C may be distinct from both Fe 
and Mg (e.g. Sanchez-Blazquez et al. 2006a; Clemens et al. 2006; 
Smith et al. 2009b), we prefer not to use Fe4668 here. Among 
the remaining indices, Fe4383 and Fe5015 are the most promis- 
ing, with good sensitivity to Fe but not to the light elements (see 
Tripicco & Bell 1995). 

Prior to describing the galaxy-by-galaxy grid inversion results, 
we present a visual impression of the data in the space defined by 
the Mgb5177, Fe5015 and HgF indices. Figures[T]and[2]show pro- 
jections of the model grid and the data, divided into six bins ac- 
cording to the velocity dispersion. The first bin contains the galax- 
ies with unresolved velocity dispersion, while the remaining bins 
represent roughly equal fractions of the galaxies with measured a. 
In Figure[T] the Mgb5177 vs Fe5015 plane is displayed; the loca- 
tion of a galaxy in this diagram reflects primarily its metallicity and 
a/Fe, but also to a lesser extent its age. The over-plotted grids are 
extracted from the TMBK models, with age chosen as appropriate 



for each a bin (determined by consideration of Figure[2j. From Fig- 
ure!]] the systematic increase of both Z/H and a/Fe with increasing 
a is apparent, although there is substantial intrinsic scatter within 
each bin. In particular, in each bin there are galaxies spanning the 
full range 0.0-0.5 in [a/Fe] covered by the models. The age effects 
are evident in Figure [2] which shows HgF versus Fe5015. Again, 
the models are shown in two-dimensional projection with the third 
variable, in this case [a/Fe], chosen to match the average value in- 
ferred from Figure[JJ The systematic trend of Z/H is visible as the 
tendency for the cloud of points to shift along the long axis of the 
grid with increasing velocity dispersion. A trend of increasing age 
is also seen in the averages. Within each bin there is substantial 
scatter, with galaxies spanning the age range 3-15 Gyr, clearly in 
excess of the measurement errors. In some of the panels it can be 
seen that the galaxies which are younger than average also have 
higher metallicity than average, i.e. there is an anti-correlation of 
age and metallicity at fixed velocity dispersion (Trager et al. 2000). 
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Figure 2. Index-index planes showing the age-sensitive index pair Fe5015 and HgF. Each panel shows galaxies within a given interval in velocity dispersion. 
For clarity, we show omit the ~4per cent of galaxies having S/N < 20. The model grids are from TMBK, with zero-points re-calibrated from the most 
massive galaxies (see Section l2~3V The grids show lines of constant [Z/H] = —0.33, 0.00, +0.35, +0.67 (blue) and of constant age = 3, 6, 9, 12, 15 Gyr (red). 
For this comparison, the grids are drawn for a fiducial a/Fe for each bin, as indicated in the legend. The fiducial a/Fe is estimated from visual examination of 
Figure[T] 



2.3 Grid inversion 

The formal grid inversions are performed using a non-linear \ 2 
minimisation over the six indices Mgb5177, Fe5015, Fe4383, 
Hbeta, HgF, HdF The predicted values for each index are rep- 
resented by a cubic-spline surface in the log(Age)-[Z/H]-[a/Fe] 
space, passing through all grid points and interpolating between. 
Using a spline interpolation ensures continuous derivatives at cell 
boundaries; the spline is constrained to reduce to linear extrapo- 
lation beyond the grid limits. Because the model grids are tilted 
with respect to the axes defined by the indices, the errors in recov- 
ered age and metallicity are correlated. We extract the full error 
covariance matrix from the \ surface, in order to propagate the 
correlated errors. 

Recalling that the indices were measured on (approximately) 
flux-calibrated spectra, and that no attempt was made to match 
the Lick system response curve, we expect some small system- 
atic offsets between the models and the observed indices. Because 
the galaxies are all at a common distance, the indices are always 
measured at approximately the same observed wavelength, and any 
offsets due to different flux calibrations should thus be similar for 
all galaxies. If uncorrected, the offsets would generate zero-point 
shifts in the relations we aim to study, e.g. the Z/H-vs-cr relation. 



A second-order effect is that, because the model grids are not per- 
fectly parallel and linear, especially over large ranges in age and 
metallicity, the slopes of the correlations could be affected by in- 
dex offsets. 

A common approach is to observe stars from the original Lick 
stellar library to establish corrections to the measured indices (usu- 
ally limited to a constant zero-point shift). This is observationally 
expensive, and overlooks the possibility of zero-point uncertainties 
in the model predictions, which may be comparable to the observa- 
tional offsets. Instead, we adopt a practical strategy of calibrating 
such that the most massive galaxies, on average, yield stellar pop- 
ulation parameters that are similar to those obtained in previous 
studies. Specifically we apply a rigid shift to the model grid, such 
that the median index-set for the 24 galaxies with a > 200 kms" 
corresponds to the prediction for (Age, [Z/H], [a/Fe]) = (10.8 Gyr, 
0.24, 0.28). These values are taken from the high-cr bin of Nelan et 
al. (2005), but they are characteristic of most other studies for gi- 
ant ellipticals, including carefully Lick-calibrated work. The index 
shifts, applied to the models (not the observations), are: (Fe5015: - 
0.69), (Hbeta: -0.10), (HdF: -0.01), (HgF: -0.05), (Mgb5177: - 
0.10), (Fe4383: -0.22), in the sense that the re-calibrated indices 
are in all cases smaller than those predicted by the published mod- 
els. We emphasize that only the average high-mass galaxy proper- 
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Predicted B-R Predicted B-R Predicted B-R 

Figure 3. Comparison of colours predicted from the spectroscopic parameters against the observed colours from the SOS catalogue of Mercurio et al. (2006). 
The observed colours, which are identical in all panels, are within an aperture of 4.4arcsec diameter, i.e. larger than the 2arcsec diameter spectroscopic 
aperture. In the left panel, the colours are predicted via the SSP models of Maraston (2005), using the measured age, but assuming a common metallicity of 
[Z/H]=0. 1 1 (the median value over the sample). In the centre panel, the measured ages are used, but a common age (6.5 Gyr) assumed. The right panel uses 
the measured values for both age and metallicity to predict the colours. The solid line shows equality; the spread is is quantified by the standard deviation (sd) 
and more robustly by the 68 per cent interval (iqr). 
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Figure 4. Residual of the observed colour, relative to the colour predicted 
from the spectroscopic age and total metallicity, compared to the spectro- 
scopic estimate of [a/Fe]. The dashed (blue) line shows a fit to the data 
after iterative rejection of outliers (crosses, without error bars). The most 
discrepant outlier NFPJ132802. 6-314521 lies beyond the plot limits. The 
dotted (red) line of slope —0.2 represents predictions from Coelho et al. 
(2007) for the effect of a enhancement on SSP colour. The vertical er- 
ror bars are dominated by contributions from the (anti-correlated) age and 
metallicity errors in the predicted colours. 

ties have been fixed in this way. Relative differences in their popula- 
tions are preserved, so that measurements of the age scatter at high 
mass, for instance, are meaningful. Moreover, the relative proper- 
ties of high- and low-mass galaxies are of course maintained, so 



that the slopes of the stellar population scaling relations are un- 
changed, to first ordei0. The calibration scheme involves only six 
degrees of freedom, compared to ~ 200 x 6 independent data. 

In Table Q] we present the derived stellar population param- 
eters from the six-index inversion, as used in analysing the pop- 
ulation trends in the following sections. Table Q] includes the fit 
results for all 232 galaxies from Paper I. For the analyses in the 
following sections, we exclude ten of these galaxies which have 
low signal-to-noise spectra (S/N < 20 per A at rest-frame 4400- 
5400 A). Velocity dispersion measurements are required for the 
analysis in Section [3j removing objects with unresolved a leaves 
193 galaxies for use in the fits. The typical formal errors are 14 per 
cent in age, 0.05 dex in Z/H and 0.04 dex in a/Fe, for galaxies 
with a — 50 — 100 km s -1 , and half of these values for galax- 
ies with a > 150 km s -1 . As a function of signal-to-noise ra- 
tio, the errors are well represented by erz/H] = 2.5(S/N)~ 1 , 
eiog(Agc) = 2.7(S/N)~ 1 , and £[Q/Fc] = 2.0(5'/iV)- 1 . 

To assess the systematic uncertainties associated with our 
choice of indices, we have investigated how the results vary when 
different sets of three indices are used to invert the model grids. In 
all cases the index set includes Mgb5177, one Fe-sensitive index 
(Fe5015 or Fe4383) and one Balmer index (Hbeta, HgF or HdF), 
giving six possible triplets. Relative to our default six-index inver- 
sion, these triplets yield zero points differing by ^-0.10 dex in age, 
< 0.05 dex in Z/H, and < 0.03 dex in a/Fe. The differences be- 
tween the results do not depend significantly on a or luminosity, 
hence the slopes of the scaling relations are quite robust. 

2 As a robustness test, we have also estimated the stellar population pa- 
rameters after applying equivalent calibrations based on the Nelan at al. av- 
erages in bins at lower masses. Although this introduces overall shifts in the 
derived parameters (e.g. 5-10 per cent in the age zero-point), we find that 
the effect of this change in calibration is not systematically correlated with 
luminosity or velocity dispersion. This test confirms that the coefficients of 
the scaling relations derived in Section[3]are stable (to within < 0.01) with 
respect to the details of the calibration scheme. 
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Table 1. Model inversion results. The first three columns give the galaxy identifier, log R-band luminosity, log velocity dispersion (cr in km s — 1 ). The following 
three columns are the log age (in Gyr), metallicity ([Z/H]) and ct-element enhancement ([a/Fe]), obtained by inversion of Thomas et al. (2003, 2004) models 
using the Mgb5177, Fe4383, Fe5015, HdF, HgF and Hbeta indices. The final three columns give the error correlation coefficients C, from which the covariance 
matrix on the stellar population parameters can be reconstructed. 



Galaxy ID 


log is. 




lof CT 
lu 6 u 




rz/Hi 

t" L L l 






[a/Fe] 


C(age,Z) 




C(Z,a) 


NFPJ132328.9-314242 


9.98 


2 


.067 ±0.019 


0.753 ±0.048 


±0.234 ±0. 


.049 


±0. 


.229 ±0.028 


-0. 


.85 


±0.32 


-0.18 


NFPJ132330.8-314935 


10.24 


2 


.014 ±0.012 


0.557 ±0.030 


±0.187 ±0. 


.030 


±0. 


.144 ±0.017 


-0. 


.85 


±0.29 


-0.16 


NFPJ132335.5-315201 


9.73 


2. 


.005 ± 0.029 


0.646 ± 0.067 


±0.261 ±0. 


.068 


-0. 


.016 ± 0.039 


-0. 


.86 


±0.16 


±0.01 


NFPJ132337.1-315047 


9.96 


2 


.129 ± 0.017 


1.041 ±0.035 


-0.080 ±0. 


.036 


±0. 


.258 ± 0.041 


-0. 


.76 


±0.42 


-0.31 


NFPJ132345.0-314230 


9.91 


1. 


.843 ± 0.039 


0.469 ± 0.060 


±0.009 ±0. 


.061 


±0. 


.061 ± 0.037 


-0. 


.93 


±0.42 


-0.33 


NFPJ132348.3-314953 


9.53 


1. 


.810 ±0.038 


0.702 ±0.062 


±0.156 ±0. 


.074 


±0. 


.100 ±0.043 


-0. 


.85 


±0.26 


-0.17 


NFPJ132355.5-313847 


9.29 


1 


.861 ± 0.046 


0.867 ±0.102 


±0.115 ±0. 


.063 


±0. 


.150 ±0.048 


-0. 


.74 


-0.03 


±0.04 


NFPJ132406.9-314449 


9.87 


2 


.041 ± 0.019 


0.985 ± 0.039 


-0.072 ±0. 


.033 


±0. 


.231 ± 0.036 


-0. 


.76 


±0.26 


-0.16 


NFPJ132412.7-314658 


9.85 


1 


.718 ± 0.074 


0.454 ± 0.072 


±0.335 ±0. 


.060 


±0. 


111 ± 0.034 


-0. 


.87 


±0.19 


-0.06 


NFPJ132418.2-3 14229 


10.31 


2 


.297 ± 0.003 


0.901 ±0.028 


±0.283 ±0. 


.021 


±0. 


.239 ± 0.011 


-0. 


.86 


±0.22 


-0.11 


NFPJ132423.0-313631 


10.05 


2 


.128 ± 0.011 


0.779 ±0.040 


±0.168 ±0. 


.035 


±0. 


.167 ±0.019 


-0. 


.87 


±0.31 


-0.21 


NFPJ132425.9-314117 


10.36 


2 


.156 ±0.007 


0.825 ± 0.048 


±0.179 ±0. 


.031 


±0. 


.161 ± 0.017 


-0. 


.83 


±0.10 


±0.00 


NFPJ132426.5-315153 


10.16 


2 


.302 ± 0.006 


0.874 ± 0.037 


±0.265 ±0. 


.029 


±0. 


.254 ±0.014 


-0. 


.87 


±0.14 


-0.07 


NFPJ 13243 1.0-3 15625 


9.56 






0.438 ± 0.082 


-0.003 ±0. 


.097 


±0. 


.022 ± 0.059 


-0. 


.94 


±0.36 


-0.23 


NFPJ132440.7-314444 


9.35 






0.595 ± 0.085 


-0.054 ±0. 


.065 


±0. 


.101 ± 0.057 


-0. 


.88 


±0.29 


-0.19 


NFPJ132507.5-314625 


9.34 






0.729 ±0.086 


-0.270 ±0. 


.055 


±0. 


.178 ± 0.078 


-0. 


.87 


±0.43 


-0.27 


NFPJ132514.9-314154 


9.85 


2. 


.319 ±0.006 


1.192 ±0.032 


-0.146 ±0. 


.022 


±0. 


.295 ± 0.024 


-0. 


.88 


±0.56 


-0.55 


NFPJ132517.1-314312 


9.51 


1. 


.880 ±0.030 


0.602 ± 0.050 


±0.189 ±0. 


.050 


±0. 


.148 ±0.030 


-0. 


.86 


±0.26 


-0.18 


NFPJ 132520.0-3 13604 


9.32 






0.307 ± 0.045 


±0.064 ±0. 


.065 


±0. 


.044 ± 0.054 


-0. 


.92 


±0.41 


-0.30 


NFPJ132526.4-314625 


9.89 


2. 


.315 ± 0.007 


1.076 ±0.016 


±0.047 ± 0. 


.021 


±0. 


.335 ± 0.017 


-0. 


.72 


±0.20 


-0.13 



2.4 Predicted colours: a consistency test 

In this section, we consider whether the stellar populations inferred 
from spectroscopy are consistent with the broad-band colours ob- 
served for the same galaxies. We use B — R colours from the 
Shapley Optical Survey (SOS) of Mercurio et al. (2006), measured 
within a fixed aperture of 4.4arcsec diameter. To transform the 
spectroscopic age and metallicity to a predicted colour, we use the 
linear relation 

(B - i?) prod = 1.150 ± 0.329 log(Age) ± 0.297[Z/H] , (1) 

derived from a fit to Maraston (2005) models with age greater than 
1 Gyr and [Z/H]> —2. Figure[3]presents comparisons between ob- 
served and predicted colours, demonstrating the comparable roles 
played by age and metallicity in reproducing the colours (see also 
Gallazzi et al. 2006). 

The measured colours exhibit a lcr scatter of 0.09 mag (esti- 
mated directly from the 68 per cent range of the data, to reduce the 
impact of outliers). If only the measured age is used to predict the 
colour (adopting a common value for the metallicity of all galax- 
ies), the residual scatter is reduced to 0.07 mag. A similar scatter 
is obtained predicting colour from metallicity and using a common 
age for all galaxies. When the measured ages and metallicities are 
used together to predict the colours, a much tighter correlation is re- 
covered, with the residual scatter reduced to 0.04 ma£l The (anti- 
correlated) errors in the age and metallicity estimates account for 
much of this scatter, leaving only ~0.02 mag unexplained. 



3 The most extreme outlier, with B — R some 0.4 mag bluer than expected 
for its age and metallicity is NFPJ132802. 6-314521. This galaxy was iden- 
tified by Miller (2005) as a powerful radio source with probable head-tail 
morphology, and detected as an X-ray source by Tajer et al. (2005). Exam- 
ination of the SOS data confirms that its central colour is indeed unusually 
blue, even compared to its large-aperture colours. 



Our predicted colours do not account for a/Fe variations, since 
Maraston's models are for solar abundance ratios. Figure [4] shows 
that the residuals from the predicted vs observed colour relation are 
mildly anti-correlated with measured a/Fe. The slope, after reject- 
ing afew outliers, is S(B-R)/S[a/Fe] = — 0.14±0. 02, indicating 
that a-enhanced populations are slightly bluer than solar-mixture 
populations. This result is qualitatively consistent with theoretical 
modelling by Coelho et al. (2007), who predict a slightly larger 
colour change 5(B — R) /S[a/Fe] ~ —0.2 mag, when Z/H is held 
constanQ. Allowing for the observed dependence on a/Fe, and for 
measurement errors, the colours show a residual intrinsic scatter 
of only 0.01 mag, compared to predictions from the spectroscopic 
parameters. 



2.5 Constraints on composite stellar populations 

We have so far been comparing observed galaxy spectra to ide- 
alised simple, i.e. single burst, stellar population (SSP) models. In 
a hierarchical galaxy formation model, real star-formation histories 
(SFHs) are expected to be extended, especially before incorpora- 
tion into a cluster, and perhaps to have multiple bursts as new gas 
becomes available from mergers. In this case, we can still fit for an 
SSP-equivalent population, which best matches the observed pa- 
rameters, but we need to be very careful regarding its interpretation. 
In general, the SSP-equivalent age is not a good approximation to 
the luminosity-weighted age. Rather, because the Balmer lines re- 
spond non-linearly with age, the SSP-equivalent age is even more 
heavily weighted to the youngest stars present. An excellent discus- 
sion of these issues, with extensive simulations of two-burst SFHs, 



4 Coelho et al. tabulated colour responses to a/Fe variations at fixed Fe/H. 
When these are transformed to responses at fixed Z/H, the increase in a/Fe 
is accompanied by a reduction in Fe/H; both effects generate bluer colours. 
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Figure 5. Comparison of age estimates using different Balmer indices. The 
upper panel shows the estimate based on HdF versus that using Hbeta, while 
the lower left plots age difference against log a. Small grey points are the 
individual data; large black points represent medians in a over the same bins 
as in Figure[T] In the lower panel, galaxies with unresolved velocity disper- 
sion have placed to the left of the dotted line. If the younger ages at low 
mass were due to "frosting" by a low-mass-fraction young sub-population, 
the data would drift to lower HdF ages at low a. In fact a slight shift in the 
opposite direction is seen. 



is given by Serra & Trager (2006); the problem was discussed pre- 
viously by Leonardi & Rose (1996) and Trager et al. (2000), among 
others. In this section, we briefly assess the evidence for composite 
stellar populations in the Shapley supercluster sample. Allanson et 
al. (2009) provide a more general analysis of our data using com- 
plex SFH models. 

Using line indices alone, it is very difficult to distinguish be- 
tween complex histories which have no star-formation in the past 
~1 Gyr, because changes in the spectrum are fairly linear in this 
regime. As a result, two widely separated star-formation events can 
produce spectra that are indistinguishable from a single burst with 
intermediate age. For such cases, a single SSP-equivalent age still 
provides a convenient moment of the age distribution, that can be 
readily compared to galaxy formation models (e.g. Smith et al. 
2009a; Trager & Somerville 2009). By contrast, sub-populations 
younger than 1 Gyr show spectral signatures of A-stars, which can- 
not be mimicked by intermediate-age SSPs. 

As a test for "frosting", i.e. young sub-populations represent- 
ing low fractions of the total mass, we can compare the SSP- 
equivalent ages determined using different age indicators. The lu- 
minosity contribution from any younger sub-population would be 



stronger at KS than at H/3, for example, due to the bluer colour 
of their main-sequence turn-off stars. As a result the H(5-derived 
ages would be younger than the H/3-derived ages. In Figure [5] we 
compare the ages derived from an inversion using Hbeta-Fe5015- 
Mgb5177 to those from HdF-Fe5015-Mgb5177. The figure shows 
that the HdF ages are in fact older than the ages based on Hbeta, 
by ~0.07 dex, the opposite of the offset expected from frosting. 
Moreover the offset is larger for low— a galaxies than for the most 
massive objects, so an increased incidence of frosting is apparently 
not the cause of the overall trend towards younger age at low mass. 
(This is also seen in the fits in Paper I, where the slope of the age-a 
relation is steeper for the fits based on Hbeta.) 

An alternative indicator for young sub-populations is the Rose 
(1984) Call index, based on the residual flux at the core of the 
Ca K and H lines, the latter being coincident with the higher-order 
Balmer line He. The index has a very non-linear behaviour, being 
approximately constant for all stars cooler than spectral type F2, 
but changing steeply for hotter stars. As a result, in composite pop- 
ulations, Ca II is very sensitive to the presence stars of < 1 Gyr 
age even if they contribute only a few percent of the mass, but lit- 
tle affected by variations in age beyond this threshold. In combi- 
nation with other Balmer indices, the Call index can distinguish 
intermediate-age approximate SSPs from old-plus-young compos- 
ite populations (Leonardi & Rose 1996). 

In Figure [6] we show Call versus the age derived from our 
default fits. For clarity we show only the 149 galaxies (64 per 
cent of the sample) with errors less than 0.15 in Call. For com- 
parison we show predictions for Call for SSPs and for popula- 
tions with a frosting of young stars, derived by combining SSP 
models from Vazdekis et al. (2009). The cases shown are for a 
~13 Gyr base population, with a secondary burst ages of 0.5, 0.7 
and l.OGyr. The mass fraction in the burst is 0.5, 1.0 and 5.0per 
cent. Both populations are assigned solar Fe/H (metallicity effects 
on Call are small enough to neglect for this test). For the hori- 
zontal axis, we compute the SSP-equivalent ages for the compos- 
ite two-burst spectrum, by inverting the TMBK models using the 
same six indices that were used for the real data. (The mixing of 
different model sets is clearly undesirable, but cannot be avoided 
for this test; the results are probably fairly generic.) These simula- 
tions predict easily-observable Call offsets for galaxies with SSP- 
equivalent ages < 5 Gyr, if their ages are driven by minor bursts 
within the last < 1 Gyr. It is seen that a few galaxies indeed lie in 
the region occupied by the frosting models, and inconsistent with 
the cool-star Call value of ~ 1.1. However, the overall distribution 
of points indicates that frosting is not responsible for the young de- 
rived ages on averag^B In particular, for 22 galaxies with inversion 
ages less than 4Gyr, the mean Call is 1.132±0.026, indistinguish- 
able from the value 1.135±0.005 for 40 galaxies with inversion 
ages above 10 Gyr (the median values are 1.128 and 1.134 respec- 
tively). 

In conclusion, we find no evidence that SSP ages are on the 
whole generated by frosting by small mass-fractions of very recent 
star-formation, although we cannot exclude this possibility for a 
small number of individual sample galaxies. The same conclusion 
holds for any other contribution from A-type stars, such as blue 
Horizontal Branch stars from an old, metal-poor sub-population 



° If the galaxies with larger Ca II errors are included, the fraction of frost- 
ing candidates remains small. This is expected, since frosted populations 
have more flux in the blue, and are thus over-represented in the low-error 
subset. 
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Figure 6. The Rose Ca II index for galaxies with errors smaller than 0.15, 
compared to the SSP-equivalent age. The blue grid indicates the expected 
behaviour for frosting by secondary bursts of age 1 .0, 0.7 and 0.5 Gyr (solid 
lines), and mass-fractions 5, 1 and 0.5 per cent (dotted lines), with the re- 
maining mass in a 1 3 Gyr base population. The red track shows predic- 
tions for 2-12 Gyr SSPs, demonstrating the stability of Ca II for intermedi- 
ate ages. Although a few galaxies fall in the region of the frosted models, 
on average the SSP-equivalent ages of "young" galaxies are not driven by 
secondary bursts in the past Gyr. 

(e.g. Maraston & Thomas 2000). Allanson et al. (2009) show that 
frosting models (or exponential decay models, which similarly al- 
low star-formation within the past Gyr) are inconsistent with the 
dynamical mass-to-light ratios implied by the Fundamental Plane. 

2.6 Corrections for stellar population gradients 

The 2 arcsec diameter fibres enclose a larger fraction of the total 
light for small galaxies than for large ones. For galaxies in our sam- 
ple, the effective (i.e. half-light) radii (R c ) range from 1 arcsec to 
10 arcsec (Gargiulo et al. 2009). Since red galaxies in general ex- 
hibit internal metallicity gradients (e.g. Davies, Sadler & Peletier 
1993), with higher metallicity in the central regions, a correction 
is needed to translate the central measurements to estimates of the 
metallicity at a common radius in terms of R c . Neglecting this cor- 
rection could lead to systematic trends that can mimic correlations 
with galaxy mass. 

To estimate the required corrections, we will assume metal- 
licity gradients of order A[Z/H]/A log R = -0.14 (Rawle et al. 
2008b), for all galaxies. Weighted by an R 1 ^ luminosity profile, 
and integrated within circular apertures, this yields a correction of 
~0. 1 dex in Z/H between a measurement within 0.1 R c and what 
would be observed within 1 R c (Figure|7]l. 

Accounting for metallicity gradients is especially important 
when we attempt to disentangle the effects of luminosity and ve- 
locity dispersion on the stellar populations. The Fundamental Plane 
implies that at fixed velocity dispersion, more luminous galaxies 
are larger. In turn this means the offset between metallicity at R e 
and metallicity within the fibre is larger for more luminous galax- 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 



Figure 7. Aperture effects on fibre measurements of [Z/H]. The curves 
shows the difference between [Z/H] measured within a circular fibre of ra- 
dius and that measured at the effective radius, R c , as a function of the 
ratio between these radii. The solid curve is for an i? 1 / 4 luminosity profile, 
and the dashed curve is for an exponential profile. We assume a differential 
metallicity gradient of 0.14 dex per decade in radius. For our Shapley sam- 
ple, the effective radii are in the range 1.0-10.0 arcsec, and i?fj b =l arcsec. 
Over this range, indicated by the vertical dotted lines, a linear approxima- 
tion to the aperture correction is justified. 

ies. Conversely, at fixed luminosity, low-er galaxies are larger than 
high-a galaxies. If uncorrected, metallicity gradients would thus 
artificially boost the correlation of Z/H with luminosity at fixed a 
and suppress its correlation with a at fixed luminosity. 

Effective radii have not been measured for all of our sample 
galaxies, but for the subset of 141 objects studied by Gargiulo et al. 
we find 

with alcr scatter of 0. 13 dex. Applying this relation (simply a state- 
ment of the Fundamental Plane) to calculate a predicted effective 
radius R^! Tcd for all galaxies in our sample, we can correct the 
metallicity to R^ led via 

A[Z/H] = -0.1 log(i?r d /^flb) • 

The corrections to Z/H range from —0.09 dex to +0.06 dex, with 
median -0.02 de)H The error in the the correction, associated 
with scatter around the Fundamental Plane, is 0. 1 x 0. 1 3 ~ 0.0 1 dex, 
which is much smaller than the formal error on Z/H (typically 
~ 0.05 dex), and is neglected in what follows. The strength of the 
internal gradients, i.e. the factor —0.1 above, is uncertain to perhaps 
50 per cent. In Section[3] we consider explicitly how our results for 
Z/H would change if the adopted aperture corrections are altered 
within this range. 

The evidence for internal age gradients in red galaxies is in- 
conclusive, with small positive (Rawle et al. 2008b) and negative 
(Sanchez-Blazquez et al. 2006b) gradients having been reported. 

e Note that within the linear part of the correction curve (Figure |7J, the 
range of aperture corrections depends only on the range of effective radius 
or equivalently the range of luminosity, if considering fixed it. It is indepen- 
dent of the actual fibre size, which affects only the average correction. 
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Table 2. Observed scaling relation fits, adopting the default aperture cor- 
rection model. For each parameter Z/H, age, and a/Fe, we provide fits with 
velocity dispersion cr, R-band luminosity L R , and to both parameters si- 
multaneously. In fits to the residuals from one of the relations (indicated by 
A), the fixed coefficient is given in parentheses. The tabulated scatter is the 
estimated intrinsic dispersion around the fit. The first line in each section 
gives the scatter around the null model, e.g. the total metallicity scatter of 
the sample as a whole. 



Relation 


Coeff of log cr 


Coeff of log L R 


scatter 


[Z/H] 


— 


— 


0.151 


[Z/H] - cr 


+0.307±0.043 


— 


0.133 


[Z/H] - L R 


— 


+0.216±0.028 


0.131 


A [Z/H] - cr 


+0.067±0.042 


(+0.216) 


0.130 


ATZ/H1 - Lb 


(+0.307) 


+0.064±0.()28 


0.131 


[Z/H] - <T, L R 


+0.149±0.062 


+0.143±0.042 


0.129 


log(Age) 






0.194 


log(Age) - cr 


+0.437±0.053 




0.165 


log(Age) - L R 




+0.112±0.041 


0.191 


Alog(Age) - cr 


+0.312±0.056 


(+0.112) 


0.176 


Alog(Age) - L R 


(+0.437) 


-0.106±0.035 


0.161 


log(Age) - cr, L R 


+0.698±0.075 


-0.235±0.050 


0.156 


[a/Fe] 






0.091 


[a/Fe] - cr 


+0.210±0.026 




0.076 


[a/Fe] - L R 




+0.047±0.020 


0.089 


A [a/Fe] - cr 


+0.156±0.027 


(+0.047) 


0.081 


A[a/Fe] - L R 


(+0.210) 


-0.053±0.017 


0.074 


[a/Fe] - cr, L R 


+0.341±0.035 


-0.117±0.024 


0.070 



Table 3. Effect of varying the aperture correction scheme on the coefficients 
obtained for [Z/H] as a function of cr and L R . 



Scheme 


AZ/H 
A log R 


Coeff of log cr 


Coeff of log L R 


None 


0.00 


+0.044 


+0.241 


Weaker 


-0.07 


+0.097 


+0.192 


Default 


-0.15 


+0.149 


+0.143 


Stronger 


-0.23 


+0.201 


+0.093 



Similarly, gradients in a/Fe appear to be small. No corrections are 
imposed here for either parameter. For the velocity dispersions, the 
usual correction is adopted, with 

, ^ prcd ) nrvll RF od 

lo § (n \ = -0.04 log — — 

(e.g. J0rgensen, Franx & Kjcergaard 1995). The corrections range 
from —0.04 dex to +0.02 dex in a. 



3 OBSERVED SCALING RELATIONS OF THE STELLAR 
POPULATIONS 

In this section we examine the dependence of age, Z/H and a/Fe 
as functions of velocity dispersion, a, and luminosity, L R . We start 
with one-parameter descriptions, e.g. the age— a relation, then test 
for residual correlations, and finally consider bivariate fits in which 
a and L R are used simultaneously to predict the stellar populations. 

We defer a discussion of the correlation among stellar popula- 
tion parameters, e.g. the age-metallicity and age-a/Fe relations, to 
a subsequent paper in this series. 



3.1 Biases and choice of the mass proxy 

Velocity dispersion has traditionally been the mass proxy used most 
widely for stellar population studies, because it is independent of 
the stellar mass-to-light ratio and reflects the depth of the gravita- 
tional potential well. However, when galaxy samples are selected 
by luminosity, the use of a as a mass proxy is biased: with decreas- 
ing a, such samples become increasingly dominated by galaxies 
with unusually high luminosity given their cr. 

The most obvious mechanism for such a bias is that old galax- 
ies are fainter than otherwise identical galaxies with younger pop- 
ulations (because stellar mass-to-light ratio increases with age) and 
are hence lost at low velocity dispersion in a luminosity-selected 
sample (as discussed by Kelson et al. 2006, for example). If not 
corrected for, this generates an artificial steepening of the age— a 
relation. To a lesser extent, the same effect steepens the Z/H— a 
relation, since the stellar mass-to-light ratio also increases slightly 
with metallicity. 

A second cause of bias can be identified, which is distinct from 
the above. Consider a parameter P, which is positively correlated 
both with stellar mass (and hence luminosity), and with a. If we 
construct the P — a relation from a luminosity-selected sample, 
the observed low-rr galaxies are unrepresentative of the population 
at large, being unusually bright for their cr. In turn, they have un- 
usually high values of P, and the slope of the P — a relation will 
be artificially flattened. This is irrespective of any direct causality 
between P and L R . A bias of this type was considered by Graves 
et al. (2007) when deriving composite spectra for SDSS galaxies 
binned by velocity dispersion. They found that within a given a in- 
terval, certain metal-line indices showed a strong correlation with 
luminosity. As a result, at low a their sample was biased to objects 
with unrepresentatively high metallicity, and they introduced a cor- 
rection method to account for the bias. In the following section, we 
will instead include the apparent luminosity dependence explicitly 
in our analysis of the scaling relations. 



3.2 Correlations with luminosity and velocity dispersion 

We begin by considering qualitatively the distribution of measured 
stellar population parameters, and their correlation with either the 
R-band luminosity L R or with velocity dispersion cr. We use only 
the set of 193 galaxies which have reliable stellar population fits 
and resolved velocity dispersions. When fitting the correlations, we 
weight galaxies according to the quadrature sum of their measure- 
ment errors and the intrinsic scatter in the fit. In practice, this is 
done by iteratively re-weighting the points, updating the estimate 
of intrinsic scatter in each step. The errors in predictors (a, L R ) 
are ignored in the fits. Table [2] presents the coefficients of the fits 
described in this section. 

The correlations with luminosity are shown in the upper pan- 
els of Figure [8] The metallicity Z/H exhibits a very strong (~8a) 
correlation with L R , while age and a/Fe are more weakly corre- 
lated (< 3cr). Because the sample is luminosity-limited, the selec- 
tion cut is vertical in these diagrams, and does not explicitly bias 
the slopes of the relationships. However, because stellar mass-to- 
light ratios depend on age, and to a lesser extent on metallicity, 
lines of constant stellar mass run diagonally across the panels. In 
the lower panels of Figure(8] we show the correlations with central 
velocity dispersion. The Z/H— a relation is slightly weaker than 
the Z/H— L R relation (~7<r), while the relations for age and a/Fe 
are much steeper (~8<r) than the equivalent L R correlations. The 
trends with a are in good agreement with the estimates made in Pa- 
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per I based on the slopes of the index-<r relations (i.e. Z/H oc a , 
age oca 052 , a/Fe oca ' 23 from table 7 of Paper I), although the 
new age slope is slightly shallower. We note that the fit to the 
a/Fe—a relation seems to be flattened by the presence of high-a/Fe 
outliers at relatively low a. The "core" of the relation would favour 
a slope ~50 per cent steeper, which would be in closer agreement 
with previous work (slopes of 0.3-0.4 were derived by Nelan et al. 
2005, Thomas et al. 2005, Bernardi et al. 2006, Graves et al. 2007). 
The fairly shallow a/Fe— a relation derived in Paper I is possibly 
influenced in the same way. 

Figure [9] shows the residual correlations. In the upper panels, 
we plot the residuals from the Z/H— a, age— a, and a/Fe—a rela- 
tions, as a function of luminosity. Fairly weak residual luminosity 
trends (2 — 4a) are found for all three parameters: Z/H increases 
with luminosity at fixed a, while age decreases with luminosity at 
fixed a; a/Fe also formally decreases with luminosity, although the 
"core" of the distribution is fairly flat. The lower panels of Figure|9] 
show, conversely, the residual dependence on a, after controlling 
for luminosity. The correlation of Z/H with a is reduced to <2a 
significance after accounting for its variation with Lr. On the other 
hand, the age and a/Fe correlations with a are strengthened when 
we control for the luminosity trend. These tests confirm the indica- 
tions from Figure[8]that velocity dispersion is dominant in driving 
age and a/Fe, while Z/H depends at least partly on luminosity. 

To explore these correlations further, we next consider bivari- 
ate scaling relations, where age, Z/H and a/Fe are expressed explic- 
itly in terms of both Lr and a, i.e. P = a log a + 61og(L_f{) + c, 
with P = log(Age), [Z/H], [a/Fe]. This formulation is unbiased by 
selection in luminosity, and also by selection in velocity disper- 
sion (recall that some 15 per cent of the Shapley galaxies are "un- 
resolved", indicating that such a selection does affect our sample). 
The robustness holds so long the distribution of galaxy properties 
at given a and given Lr is not affected by additional selection cri- 
teria. 

We present the coefficients of the bivariate fits in Table [2] 
and illustrate the results graphically in Figure [10] For both age and 
a/Fe, the bivariate fits confirm a positive correlation with a and si- 
multaneously a negative correlation with luminosity. For Z/H, there 
are positive dependences on both a and L. Table [3] shows that the 
luminosity correlations are sensitive to the adopted aperture cor- 
rections. If no correction at all is allowed for the metallicity gra- 
dients, the recovered trend of Z/H would be primarily with lumi- 
nosity. An alternative way to present the dependence of parameters 
on luminosity and velocity dispersion is through "marked Faber- 
Jackson" relations (Figure [TTt in which the location in the (L, a) 
plane is colour-coded by age, Z/H, a/Fe in turn. These figures pro- 
vide an immediate visual confirmation of the different behaviour of 
Z/H (inclined contours) relative to the other parameters (contours 
aligned mainly with a). A somewhat comparable illustration, for 
age, is shown by Gallazzi et al. (2006) in their figure 15a. 

3.3 Bivariate scaling relation for the broadband colour 

Using the bivariate scalings of age and metallicity with velocity dis- 
persion and luminosity, it is now possible to predict the equivalent 
correlations for broadband colours, and compare to the observed 
trends. 

A bivariate fit to the observed B — R colours from SOS, within 
an aperture of 4.4arcsec diameter, yields a correlation primarily 
with velocity dispersion: A(B - R) ohs = (0.19 ± 0.03) A log a - 
(0.01 ± 0.01) A log L. This agrees with the findings of Bernardi et 
al. (2005) who showed that a strong colour— a relation holds after 



controlling for luminosity, but that there is no colour-luminosity 
relationship at fixed a (their Figure 4). 

We can readily explain this behaviour using our 
spectroscopically-determined age and metallicity scaling re- 
lations. Since the observed colours are within a fixed aperture, 
we use the scalings derived for the "no aperture correction" 
case. The factor-of-two mismatch between the photometry aper- 
ture and spectrograph fibre size should yield only a constant 
correction factor, and will not affect the slopes of the correla- 
tion. Using Equation [T] which was derived from the Maraston 
(2005) models, predicted colour differences can be expressed 
as A(B - H) pred = 0.33A log(Age) + 0.30A[Z/H]. The 
small effects of a/Fe on broadband colours will be neglected 
here. Inserting the L, a dependences of age and metallicity, this 
gives A(B - i?) prcd = 0.23Alogcr - 0.08AlogL from the 
age change and A(B - i?) prcd = 0.01 A log a + 0.07A log L 
from Z/H variations. We see that the age and metallicity con- 
tributions to the L dependence cancel almost precisely, leaving 
A(B - R) prc d = 0.24A log a - 0.01A log L, dominated by the 
change in a, as observed. 

This test confirms the consistency between the derived spec- 
troscopic parameters and the observed broadband colours. Where 
Section l2T4l demonstrated this on a galaxy-by-galaxy basis, here we 
have shown that the derived systematic behaviour is also reflected 
in the colours. 



4 NO CORRELATIONS WITH STELLAR MASS? 

We noted in Sec tion [3Tl that simple fits for correlations with a are 
compromised for samples selected by luminosity, and that although 
introducing luminosity as a second "predictor" leads to unbiased 
fits, these are difficult to interpret because luminosity itself depends 
on the stellar populations. A physically more meaningful quantity 
to use as a predictor of age, etc, is the stellar mass, M* . In this sec- 
tion, we develop a Monte Carlo modelling approach to constrain 
the relative importance of M* and a in driving the stellar popula- 
tion scaling relations. 

4.1 Modelling the scaling relations 

In principle, the stellar mass-to-light ratio M»/L can be derived 
from the spectroscopic age and metallicity estimate^. However, 
it is not straightforward to use stellar masses derived this way to 
fit the scaling relations, because the ages and stellar masses will 
have near-perfectly correlated measurement errors. A more general 
obstacle is that at low stellar mass, the sample is necessarily biased 
by the luminosity selection towards younger and brighter galaxies. 

To estimate the correlations of age and metallicity with stel- 
lar mass, we instead resort to modelling the assumed underlying 
distribution of galaxy properties such that a simulated luminosity- 
selected sample approximately reproduces the data. In practice, this 
modelling requires a priori unknowns such as the distribution func- 
tion of stellar masses and the dependence of velocity dispersion 

7 For samples of galaxies having a wide intrinsic range in M, /L, optical 
colour and near-infrared luminosity is widely used to estimate stellar mass 
(e.g. Bell et al. 2003). This method successfully discriminates between low- 
mass star-forming galaxies and more massive but passive objects. By con- 
trast for old SSP models, i.e. within the passive galaxy population, optical 
colour is a very poor predictor of M* /Ljf , and hence this approach does 
not provide an improved estimate of M* . 
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Figure 8. Stellar parameters versus two alternative mass proxies. The upper panels show the trends with luminosity, i.e. the Z/H— Lr, age— Lr and a/Fe— Lr 
relations. The sample selection limit is indicated by the vertical dashed line. For the Z/H panel, lines of constant stellar mass (M* = 10 9 ' 5 , 10 10 ' 5 , 10 11 ' 5 Mq) 
are shown by dotted lines, assuming an age of lOGyr. In the age panel, solar metallicity is assumed for the stellar mass lines. In the lower panels, we show 
the equivalent trends with velocity dispersion a. The horizontal axis scales in the lower panels have been matched to the upper panels using the relationship 
between Lr and log a for this sample. The significance of the fitted slope (in units of the standard error) is indicated in parentheses in the title-bar. For the 
coefficients of the fits, see Table|2] 



Table 4. Results for Monte Carlo modelling of the underlying scaling relations. The input model assumes that the stellar populations depend intrinsically only 
on velocity dispersion (<r), and are not systematically related to stellar mass (M„). The recovered correlations with a and luminosity (Lr) are comparable to 
the observed relations, when correlated residuals (the Z-plane) are allowed for. Solutions (1), (2) and (3) are as labelled in Figure [i"2l 



Input model <r, M* scalings 


Z/H 


oc a+°- 35 M°- 00 


Age 


oc 


CT +0.40 M 0.00 


a/Fe 


oc a+ 20 M,°- m 


no M* dependence 




Recovered a, L r scalings 


Z/H 
Z/H 
Z/H 


oca+ 0A1 L R - 05 
oc a+ 0M L+ - 21 
oc a+°- 17 L+ - 14 


Age 
Age 
Age 


oc 
oc 
oc 


^+0.78^-0.31 

CT +0.7l£-0.25 

_+0.74 r -0.28 
a ^R 


a/Fe 
a/Fe 
a/Fe 


oc a+°- 20 L R om 
oc a+ 34 L R 011 
oc a+°^L R - 10 


uncorrected age-Z residuals 
perfect anti-correlation 
moderate anti-correlation 


(1) 

(3) 
(2) 


Observed cr, Lr scalings 


Z/H 


oc CT+O-iSL+O- 14 


Age 


oc 


^+0.70^,-0.24 


a/Fe 


oc CT +0-34 L -0 l2 


(default aperture correction) 





on stellar mass. Thus the models must be tuned to reproduce both 
the luminosity distribution and the observed Faber-Jackson rela- 
tion of the galaxy sample, as well as the observed scalings of age 
and metallicity with luminosity and velocity dispersion. 

To construct a simulated galaxy sample, we first draw stel- 
lar masses from a Schechter function, with characteristic mass 
M* = 1O 1O ' 7 M0, and low-mass slope qmf = —0.5, and assign 
velocity dispersions with mean value following log a = 1.85 + 
0.50 log w iti* M ; » an d a gaussian scatter of 0.125 in log a around 
the mean. The numeric coefficients here are chosen (but not rig- 



orously fitted) to match the luminosity distribution and the ob- 
served Faber-Jackson relation of our galaxy sample. The adopted 
stellar mass function is in fact similar to that derived for early- 
type (concentration-selected) galaxies by Bell et al. (2003), where 
Mt = 10 10 - 6 M Q and a M F = -0.6. We attempt to reproduce 
the observed stellar population scaling relations P = P a a L b 
in terms of underlying models of the form P = Pq a a M\ , 
where the P are age, Z/H and a/Fe. For given choices of model 
parameters (Pg,a' P ,bp), we assign ages and metallicities, with a 
fixed scatter s p around the model. We allow the deviations from 
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Figure 9. Residual correlations with the mass proxies. In the upper row, we show the correlations of age, Z/H and a/Fe with luminosity, after removing the 
trends with velocity dispersion. In the lower panel, the trends with <r are shown after removing the luminosity trends. Note the different behaviour of Z/H 
compared to the other two parameters. The significance of the fitted slope (in units of the standard error) is indicated in parentheses in the title-bar. For the 
coefficients of the fits, see Table[2] 



[Z/H] 



log(Age) 



[a/Fe] 




0.0 0.2 0.4 
coefficient of log(a) 



0.0 0.5 1.0 

coefficient of log(o) 



0.0 0.2 0.4 

coefficient of log(a) 



Figure 10. Bivariate model coefficients. The shaded ellipses show the 1, 2, 3 a intervals on the coefficients of fits log(Age)=ao + ai log a + (22 log Lr, and 
similarly for [Z/H] and [a/Fe]. For the metallicity panel, triangles indicate the fits obtained adopting 50 per cent stronger (red, upward-pointing) and 50 per 
cent weaker (blue, downward-pointing) aperture corrections for metallicity gradients (see Table[5J. The open circle shows the fit when no aperture correction 
is applied at all. Because the predictors a and Lr are correlated, the errors on their coefficients are anti-correlated. The unshaded grey ellipses indicate the 
correlations obtained from fitting the SDSS data from Graves et al. (2009), as discussed in Section[5] 
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[Z/H] log(Age) [a/Fe] 




log(o) log(o) log(o) 

Figure 11. Marked Faber-Jackson relations, showing the distribution of SSP parameters in the space of luminosity and velocity dispersion. The colours are 
assigned according to rank in each parameter shown. The contours show lines of equal parameter rank, with a heavy smoothing applied. 



[Z/H] log(Age) [a/Fe] 




-0.2 0.0 0.2 0.4 0.6 0.0 0.5 1.0 0.0 0.2 0.4 0.6 



coefficient of log(o) coefficient of log(a) coefficient of log(a) 

Figure 12. The observed correlations (yellow ellipses, as in Figure fTol compared to the outputs from simulations that assume no intrinsic dependence on stellar 
mass. The velocity dispersion dependence of the input models is indicated by the blue square on the horizontal axis. The correlations with a and L, obtained 
by fitting the simulated dataset, are shown by the numbered red circles, for three cases described in the text: 1 for uncorrelated age/metallicity scatter; 2 for 
moderately anti-correlated scatter; 3 for perfectly anti-correlated scatter. The figure demonstrates that the observed luminosity dependences can be recovered 
without any intrinsic correlation with stellar mass, provided that the age/metallicity scatter is anti-correlated. 



the mean scaling relation for age to be anti-correlated with those 
from the Z/H scaling relation, as proposed by Trager et al. (2000). 
The a/Fe deviation is assumed to be correlated with the age de- 
viation. We consider cases of uncorrelated, moderately correlated 
and perfectly correlated deviates. Using the Maraston (2005) SSP 
models, we next determine the stellar mass-to-light ratio for each 
simulated galaxy, based on its age and Z/H, and from this com- 
pute the luminosity Lp. (Abundance ratio effects on the mass-to- 
light ratio are neglected here.) Finally, we impose selection criteria, 
logfP/Z/©) > 9.2 and logo - > 1.55 to reproduce the approximate 
limits of the observed galaxy sample. 

The P = P a a L b relations can now be fit for the selected 
simulated galaxies, and compared to the observed correlations. Ide- 
ally, the model would be optimized to obtain confidence inter- 
vals for the parameters which describe the underlying correlations 
(i.e. Pq, a' P , b' P , sp, for each parameter P), marginalizing over the 
"nuisance parameters" (i.e. those describing the stellar mass dis- 
tribution function, the M* — a relation, and the correlation struc- 
ture of the deviations from the mean scalings). In practice how- 
ever, the high dimensionality precludes a full exploration of the 



solutions, at least within this paper. Here, we simply test whether 
suitable solutions exist for underlying relations having no intrin- 
sic dependence on stellar mass, i.e. all stellar population parame- 
ters are determined only by their velocity dispersion. Table [4] and 
Figure [T2] summarize the results obtained for input correlations of 
Z/H oc a+°' 35 ,Age oc a+ ' 40 , and a/Fe oc <7+ - 20 . As we 
will show, these input slopes are able to reproduce all of the ob- 
served trends with a and L, without invoking intrinsic stellar mass 
dependence. The required slopes were determined using a three- 
parameter grid search, with the parameters of the stellar mass dis- 
tribution and the M, — a relation held fixed, and adopting the case 
with "moderate" anti-correlation of age and metallicity residuals 
(see below). The scatter in parameters around the input model is 
constrained in all cases to reproduce the age, metallicity and a/Fe 
scatter around the observed scaling relations. 
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[Z/H] vs luminosity (8.9) Age vs luminosity (3.2) [a/Fe] vs luminosity (2.2) 




log a log a log a 

Figure 13. Stellar parameters versus two alternative mass proxies, as derived from our Monte Carlo model in which there is no intrinsic dependence on stellar 
mass. The figure shows a single realisation for ~200 galaxies selected from the model to match the observed sample definition. The grey points indicate 
galaxies from the model that were excluded by the luminosity and a limits. The overall pattern of results should be directly compared with the observed data, 
in Figure[8] 



4.2 Results 

In the case where the age and metallicity deviates are uncorrelated 
(solution 1 in Figure [72l and Table [4}, we find that although age 
does not depend explicitly on stellar mass, the recovered fit has age 
increasing very steeply with a but also decreasing with luminosity. 
This is primarily an effect of the variation in M* / L: younger galax- 
ies are brighter. For Z/H and a/Fe the recovered scaling relations 
are much closer to the input relations, because M* / L depends only 
weakly on Z/H and not at all (in our modelling) on a/Fe. Thus, if we 
do not allow for correlated deviations, our models with no intrinsic 
stellar mass dependence cannot match the observed trends in these 
parameter^. Indeed, the recovered Z/H scaling would yield a lumi- 
nosity slope opposite to that observed, because higher-metallicity 
galaxies have lower luminosities, other factors being equal. 

Next we impose an anti-correlation between the age and 
metallicity deviations at fixed mass, such that galaxies that are 
younger than average are also more metal rich than average. This 
is consistent with the "Z-plane" in our own data (Smith et al. 2008) 
and with the conclusions reached in other work (e.g. Trager et al. 
2000; Smith et al. 2009a). We also allow a positive correlation of 

8 We estimate that to match the Z/H trends in the absence of correlated 
deviation, we would require approximately Z/Hoc a +0,1 M^~ ' 2 . 



the a/Fe and age deviates. In the first instance we assume perfect 
anti-correlation among the residuals. In this case, the age trends 
are only slightly altered from the uncorrelated-deviation solution. 
However, the recovered metallicity scalings are dramatically differ- 
ent: instead of the input pure a dependence, the observed scaling 
is almost purely with luminosity, in the sense that more luminous 
galaxies have higher Z/H. The explanation is that more luminous 
galaxies at given a and M, are younger, and the Z-plane now im- 
plies that they are also more metal rich. This indirect brightening 
effect is much larger than the direct dimming effect of higher Z/H. 
A similar effect holds for a/Fe. The recovered scaling relations are 
now in much better agreement with the observed correlations, al- 
though for Z/H they "overshoot", producing slightly too strong a 
luminosity trend. Adjusting the input a scaling does not much im- 
prove the match. 

Finally, noting that the observed trends lie intermediate be- 
tween the results for uncorrelated and perfectly anti-correlated age- 
metallicity deviates, we assign stellar populations to the simulated 
galaxies using moderately anti-correlated scatter (p = —0.75). 
With this prescription, we are able to match the observed scal- 
ing relations for all three stellar population parameters. Figure [T3l 
shows the simulated one-parameter scaling relations (i.e. equiva- 
lent to Figure [§}, for a realisation of our model, with moderate 
anti-correlation of the age/metallicity deviates. Exactly as in the ob- 
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served galaxy sample, we recover metallicity scaling tightly both 
with luminosity and with a, while age and a/Fe follow a much 
more closely than luminosity. We emphasise that these predictions 
are derived from an underlying model in which stellar populations 
depend only on a, and are independent of stellar mass. This re- 
sult is robust to changes in the strength of the assumed metallicity 
gradients: e.g. if stronger aperture corrections are preferred, the de- 
gree of residual anti-correlation required is reduced. Note also that 
the scalings with velocity dispersion in the Monte Carlo models 
(e.g. Age tx <t +0 ' 40 ) are not substantially different from the single- 
parameter fits of Section [X2l (where Age oc a + 0A4±0 - 05 y j n princi- 
ple, the modelling here explicitly takes into account the biases de- 
scribed in Section l3TI The similarity to the one-parameter scalings 
confirms that selection biases do not strongly influence the simpler 
fitting method (as found also by Allanson et al. 2009). 

With the simulations in hand, we can assess the apparently 
simpler approach of fitting scaling relations using stellar masses 
estimated directly from the measured ages and metallicities. As 
noted above, this is complicated by the very strong correlation be- 
tween age and spectroscopic stellar mass, and by the luminosity 
selection limit, which translates into a sharp diagonal cut in the 
age-i\f» plane. Figure [14] shows the age-M» relation for the ob- 
served sample, and for a realisation of our simulation, in which 
there is no explicit dependence on stellar mass. In both cases a 
strong correlation is observed, but this primarily reflects the corre- 
lation with velocity dispersion, and bias from the luminosity limit. 
Fitting simultaneously using a and A/*, the observed data yield 
Ageoc CT +°- 2 6±°- 1 °M+°' 13±0 - 05 . Hence naively we would infer 
a significant correlation of age with stellar mass, at fixed velocity 
dispersion. However, a similar correlation is recovered for our sim- 
ulated samples, which have no intrinsic M* dependence. Specifi- 
cally, we obtain Ageoc a + - 24 M+ - 09 w hich is indistinguishable 
from the measured trends. In both the simulated and the observed 
samples, the apparent dependence on M* (at fixed velocity disper- 
sion) arises from error correlations and from the luminosity selec- 
tion bias. 

We conclude that for reasonable assumptions about the galaxy 
population (mass distribution, M» — a relation, Z-plane anti- 
correlation of age and metallicity residuals at fixed a), we can 
reproduce the observed scaling relations from a model in which 
stellar populations depend only on velocity dispersion, and not on 
stellar mass. A residual correlation with M* can not be excluded, 
since the assumptions entering into the model are uncertain to an 
extent, and a full exploration of the parameter space has not yet 
been performed. Rather, our conclusion is that the data do not at 
present require any stellar mass dependence, in contrast to recent 
claims, which we discuss in the following section. 



5 DISCUSSION 

During preparation of this paper, Graves et al. (2009, hereafter 
GFS09) published an analysis of composite spectra from the SDSS 
which addresses some of the issues that are discussed here. At an 
observational level, the results of GFS09 agree closely with the fits 
we obtain from our sample of galaxies in Shapley. However, our 
interpretation of the results, and especially our conclusion that no 
residual dependence on A/* is required, differ from that of GFS09, 
as we discuss in this section. 

GFS09 constructed composite spectra for 16000 SDSS galax- 
ies with 0.04 < 2 < 0.08, binned by a, L and optical colour. The 
high-S/N stacked spectra for each bin were used to derive SSP- 
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Figure 14. The age versus stellar mass relation in the observed sample 
(above) compared to a representative realisation from the simulation (be- 
low), in which age oc <j+ 4 but there is no intrinsic dependence on Af*. 
For the observations, stellar masses have been estimated using the ob- 
served ages and metallicities. The dashed line in each panel indicates the 
10 9 ' 5 Mq luminosity limit, assuming a fixed metallicity. For the simulated 
data, the figure assumes realistic measurement errors in the age and metal- 
licity. The stellar masses shown are those which would be computed from 
the spectra by an observer, and hence include the effects of these measure- 
ment errors. Grey points indicate simulated galaxies excluded by the lumi- 
nosity and <r limits. 



equivalent parameters, via comparison to models from Schiavon 
(2007), which were then used to explore how the stellar popula- 
tions are correlated with luminosity and velocity dispersion. Their 
galaxy sample was selected to have concentrated 7- 1//4 -like lumi- 
nosity profiles, as well as to have no detected emission lines, con- 
trasting with the absence of explicit morphological cut in our sam- 
ple. The total range of galaxy mass is narrower in GFS09 than here, 
with loga > 1.86 and log(L/L0) > 9.4. Moreover the sample 
is drawn from the general galaxy population, rather than being ex- 
plicitly a cluster-galaxy sample as here. GFS09 did not apply any 
correction for metallicity gradients, stating that the SDSS 3-arcsec 
fibres are "not nuclear spectra but sample a substantial fraction of 
the galaxy light". 

Using the stellar population parameters for the stacked spec- 
tra, GFS09 recover positive correlations of age, Fe/H and Mg/Fe 
both with velocity dispersion and with luminosity. In agree- 
ment with our results, the correlations for age and Mg/Fe are 
much stronger when a is the predictor variable, while Fe/H is 
more strongly correlated with luminosity. GFS09 did not pro- 
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vide quantitative descriptions of these correlations. Fitting to the 
data in their tables 1 and 2, we obtain Age oc cr +0,95 L" - 25 , 
Fe/Hoc cr- 010 L +0 ' 28 and Mg/Fe oc cr +a53 L~ 012 . For a fair com- 
parison to our results for metallicity, we use the translation derived 
by Smith et al. (2009a) to convert from Fe/H and Mg/Fe in the 
Schiavon models to Z/H in the TMBK models. Fitting the result- 
ing quantity yields Z/Hoc a" 1 " ' 30 L +0 ' 21 . These results are quali- 
tatively similar to ours, although outside our formal error ellipses 
mainly due to a stronger a dependence for Mg/Fe, which affects 
both the o/Fe and Z/H comparisons. 

Although we broadly agree with GFS09 in terms of the ob- 
served correlations, we differ in our interpretation of the luminos- 
ity component of the trends. GFS09 claimed that the variations in 
M*/L due to varying age are not sufficient to generate the luminos- 
ity effects observed, and asserted that the luminosity dependence of 
the stellar populations reflects a real correlation with stellar mass, 
rather than "trivial" effects due to M»/L, In Section [4] we ex- 
plicitly examined the effects of M*/L variations through Monte 
Carlo simulations including the effects of sample selection, as well 
as intrinsic correlations between age and metallicity at fixed mass. 
We showed that, for our observational strategy, stellar populations 
which depend only on a can give rise to luminosity dependence 
consistent with that observed. It seems plausible that the similar 
trends found by GFS09 can be explained similarly, but to confirm 
this would require equivalent Monte Carlo modelling tuned to the 
different strategy used in their work. 

Moreover, we have seen that the effect of metallicity gradi- 
ents can be substantial when deriving trends with a and L simulta- 
neously, because the Fundamental Plane implies that the brightest 
galaxies at a given velocity dispersion have larger effective radius, 
and are thus subject to larger aperture corrections. The range of 
aperture corrections at fixed a, over a given interval in luminos- 
ity, is independent of fibre size and independent of distance to the 
galaxy. Thus despite GFS09 using larger fibres, and covering a sig- 
nificant range in distance, we expect that uncorrected metallicity 
gradients will generate a bias in their Fe/H scaling relations, in 
favour of a steeper slope with L, and a shallower dependence on 

(7. 

Finally, we note that our results partly disagree with those of 
Gallazzi et al. (2006), who found that the anti-correlation of age 
and luminosity at fixed a in an SDSS-based sample could not be en- 
tirely accounted for by variations in M,/L and hence some resid- 
ual dependence with Af» was required. For metallicity, however, 
Gallazzi et al. conclude that dynamical mass is more important than 
stellar mass, in accordance with our results. 



6 CONCLUSIONS 

In this paper we have used previously-reported absorption line 
strengths for passive galaxies in the Shapley Supercluster to obtain 
estimates of SSP-equivalent age, metallicity and a-element abun- 
dance ratios. Typical formal errors on the derived ages and metal- 
licities are ~ 10 per cent. The derived population parameters repro- 
duce the observed broadband colours, with a weak residual colour 
dependence on a/Fe in the sense predicted by stellar evolution 
models. Using the Rose Ca II index, we have shown that the young 
(2-5 Gyr) SSP-equivalent ages in some galaxies do not, in general, 
result from a low mass-fraction "frosting" of star-formation within 
the past Gyr. 

We investigated the variation of stellar population parameters 
as a function of both velocity dispersion and luminosity. Examin- 



ing both the residual correlations after fitting one variable, and the 
bivariate fits, we find strikingly different behaviours for the differ- 
ent parameters. Age and a/Fe are correlated primarily with a, but 
with an additional trend with luminosity such that at given a, the 
more luminous galaxies are younger and less Q-enriched. For total 
metallicity Z/H, the observed correlations appear to be primarily 
with luminosity. These observed trends, for galaxies in rich cluster 
cores, are in close agreement with the recent results of Graves et al. 
(2009) for a general sample of galaxies from SDSS. Thus it is un- 
likely that the results are strongly dependent on the environments 
that are sampled. 

The luminosity trends would naively suggest that metallicity, 
in particular, may be primarily correlated with stellar mass, rather 
than with the depth of the potential well. Such a result appears 
contrary to galactic wind models for the mass-metallicity relation. 
However, using Monte Carlo simulations, we showed that the ap- 
parent luminosity correlations can be recovered from an underlying 
model in which stellar populations depend only on velocity disper- 
sion, and not at all on stellar mass. Thus the observed luminos- 
ity trends do not imply corresponding trends of stellar population 
parameters with stellar mass. For age, the luminosity correlation 
arises because young galaxies are brighter, at fixed M* . For metal- 
licity, the observed luminosity dependence arises because metal- 
rich galaxies, at fixed mass, tend also to be younger, and hence 
brighter. The presence of metallicity gradients also boosts the ap- 
parent luminosity dependence, if aperture corrections are not ap- 
plied. 

After accounting for these effects, we conclude that the stellar 
populations of passive galaxies are fundamentally correlated with 
velocity dispersion, and hence with the depth of their gravitational 
potential well. There is no convincing evidence for additional de- 
pendence on the stellar mass, after accounting for the velocity dis- 
persion correlations. 
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